Chain conformation of ring polymers under a cylindrical nanochannel confinement 
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We investigate the chain conformation of ring polymers confined to a cylindrical nanochannel 
using both theoretical analysis and three dimensional Langevin dynamics simulations. We predict 
that the longitudinal size of a ring polymer scales with the chain length and the diameter of the 
channel in the same manner as that for linear chains based on scaling analysis and Flory-type theory. 
Moreover, Flory-type theory also gives the ratio of the longitudinal sizes for a ring polymer and a 
linear chain with identical chain length. These theoretical predictions are confirmed by numerical 
simulations. Finally, our simulation results show that this ratio first decreases and then saturates 
with increasing the chain stiffness, which has interpreted the discrepancy in experiments. Our results 
have biological significance. 

PACS numbers: 87.15.-v, 82.35.Lr, 87.15.H- 



Introduction. The ring polymer, also called cyclic poly- 
mer, is an important member of the polymer family. In 
biopolymer science, there exist circular DNA, cyclic pep- 
tides and cyclic oligosaccharides and polysaccharides 

The statistical mechanical properties of polymer are 
influenced by the ring closure. One of the fundamental 
properties for ring polymers is the chain sizes in solution. 
In a 6 solvent, based on Gaussian approximation Zimm 
and Stockmayer 0] have shown the radius of gyration of 
ideal ring polymers, R g , r ~ N v with v = 0.5, which scales 
with the chain length N in the same manner as that of 
ideal linear polymers, R g .i, and further found the ratio 
Rg.r/Rg.i = l/v2 = 0.707 for ring polymers and linear 
chains of the same chain length. Under a good solvent 
conditions, based on a renormalized two-parameter the- 
ory Douglas and Freed Q have demonstrated R g , r ~ N u 
with v = 0.6, where the scaling exponent is also the same 
as that for a linear polymer. In addition, they have also 
found the ratio i? g>r /i? Si / = 0.718, which is larger than 
that for ideal chains. The scaling exponent v = 0.6 is 
in agreement with other theoretical results 0-0] and ex- 
periments 0, but slightly larger values of R g . r /R g .i arc 
observed in simulations and experiments [l[ . 

Motivated by the fundamental relevance in polymer 
physics and many biological processes, such as DNA 
packaging inside the phage cap sid 0] , polymer transloca- 
tion through nanopores [9|, [l0( and viruses injecting their 
DNA into a host cell [ll[ , the properties of a polymer con- 
fined in a nanochannel have been the subject of extensive 
studies flil - Erjl ]. Based on the blob picture [22j|, the con- 
formational properties of a linear self-avoiding polymer 
chain confined in a slit or in a cylindrical nanochannel are 

Unlike its linear poly- 



relatively well understood (27H3C 
mer counterpart, the physics of confined ring polymers is 
still not clear. Only recently, few studies have addressed 
scmiflexible ring polymers. For semiflexible ring poly- 
mers in weak spherical confinement, the internal struc- 



ture of the chain shows bucklin g an d a conformational 
transition to a figure eight form [3l|. For the conforma- 
tional properties of a semiflexible ring polymer confined 
to different geometrical constraints, Fritsche and Heer- 



mann |32j found an important role of the geometry of 
confinement in shaping the spatial organization of poly- 
mers. 

However, even for flexible ring polymers confined to 
a cylindrical nanochannel, the behavior of the chain re- 
mains unclear. The basic questions associated with flex- 
ible ring polymers confined to a cylindrical nanochan- 
nel are the following: (a) What's the scaling behavior of 
the longitudinal size with the chain length and the chan- 
nel diameter? (b) What's the ratio of the longitudinal 
sizes for a ring polymer and the linear chain of the same 
length? The second question deals with the numerical 
prefactors of the scaling behaviors, which may change our 
understanding about polymers. On the one hand, the 
numerical prefactors is very important for understand- 
ing the biopolymer dynamics, such as the chain ejection 
out of the nanochannel where the confinement induced 
driving force is related to the numerical prcfactor. On 
the other hand, knowing numerical prefactors can help 
us deduce the chain conformation of ring polymers con- 
fined to a cylindrical nanochannel by comparison with 
the wellknown chain conformation of the confined linear 
chain. To this end, in this work we examine the proper- 
ties of ring polymers confined in nanochannel and their 
differences from the linear chains. 

Model and methods. In simulations, the polymer is 
coarse-grained as bead-spring chain of Lcnnard- Jones 
(LJ) particles with the Finite Extension Nonlinear Elas- 
tic (FENE) potential. We use a short range repulsive LJ 
potential to model the excluded volume interaction be- 
tween beads: U LJ (r) =4£[(f) 12 -(2) 6 ]+£forr < 2 1 / 6 cr 
and for r > 2 1 / 6 cr. Here, a is the diameter of a bead, 
and £ is the depth of the potential. The connectivity be- 
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tween neighboring beads is modeled by a FENE spring 
with Ufene{t) = — \kR^h\(l — r 2 /i? 2 ,), where r is the 
distance between consecutive beads, k is the spring con- 
stant and Rq is the maximum allowed separation between 
connected beads. 

The nanochannel is described as stationary particles 
within distance a from one another which interact with 
the beads by the repulsive Lennard- Jones potential. The 
nanochannel particle positions are not changed in the 
simulations. Each bead is subjected to conservative, fric- 
tional, and random forces in the Langcvin dynamics sim- 
ulation H3: mr l = -V(U LJ + U FE NE)-^i + Ff. Here 
to is the bead's mass, £ is the friction coefficient, Vj is 
the bead's velocity, and Ff is the random force which 
satisfies the fluctuation-dissipation theorem. In the sim- 
ulation, we use the LJ parameters e, o~, and to to fix the 
system energy, length and mass units, respectively. Thus, 



we have the corresponding time scale t^j = (ma /e) 



1/2 

and force scale e/a, which are of the order of ps and pN, 
respectively. The dimensionless parameters in the model 
are then chosen to be Rq = 1.5, k = 15, £ = 0.7. 

In our model, each bead corresponds to a Kuhn length 
(twice of the persistence length) of a polymer. For a 
single-stranded DNA (ssDNA), the persistence length of 
the ssDNA is sequence and solvent dependent and varies 
in a wide range [HI, S]. Here, we choose the value of 
a ~ 2.8 nm for a ssDNA containing approximately four 
nucleotide bases (35|. The average mass of a base in 
DNA is about 312 amu, leading to the bead mass to sa 
1248 amu. We set fcsT = 1.2£, and thus the interaction 
strength e is 3.39 x 10~ 21 J at actual temperature 295 K. 
This leads to a time scale of 69.2 ps and a force scale of 
1.2 pN. The Langcvin equation is then integrated in time 
by a method described by Ermak and Buckholz (36[. 

Results and discussion. One of the central properties 
of a polymer confined in a cylindrical nanochannel is the 
longitudinal size of the polymer. Consider a flexible poly- 
mer confined to a cylindrical nanochannel of diameter D 
which is less than R g , the radius of gyration of the chain. 
The Response of the polymer to confinements is primar- 
ily dictated by the relative value of D with respect to 
the chain persistence length. Depending on whether D is 
larg er (de Gennes regime 27j) or smaller (Odijk regime 
[37l | ) than the chain persistence length, different scaling 
behaviors of the longitudinal size of the chain, Rn , as a 
function of D were predicted in the pioneering theoreti- 
cal studies by de Gennes [13] and Odijk jjjll, respectively. 
In the Odijk regime, the physics is decimated not by ex- 
cluded volume but by the interplay of confinement and 
intrinsic polymer elasticity. To consider Odijk regime, 
we need to take into account the chain stiffness in the 
model. In this work, we only consider the de Gennes 
regime, namely a -C D <?C R g with a being the length of 
a segment of the chain. 

For a single polymer chain confined in a cylindrical 
nanochannel of diameter D, there are two characteristic 



lengths, D and R g . Based on the scaling analysis [271 ]. 
the length occupied by the chain in the nanochannel, Ru 
has the scaling form 



Rll Rn 



R 



D h 



(1) 



where <fr(x) ~ x m is a dimensionless scaling function. 
In the de Gennes regime, the chain becomes a one- 
dimensional chain, noting to do with the chain topology 
(ring polymers or linear chains). Therefore, for both ring 
polymers and linear chains, Rn ~ N. As a sequence, 
to = 1/v — 1 with v = 0.588 is the flory exponent in 
three dimensions j27|. According to Eq. ([T}, we have 
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where, R\\ T and Rnj are the lengths occupied by the 
chain in the nanochannel for ring polymers and linear 
chains, respectively. Although above scaling analysis 
clearly provides the physics for i?y, we still cannot ob- 



tain the ratio 



To this end, we use Flory-like theory 



for both linear chains and ring polymers. 

According to the blob picture, for a linear chian con- 
fined in an infinitely long nanochannel of diameter D, the 
chain will extend along the channel axis forming a string 
of blobs of diameter = D. On length scales smaller 
than the blob size, the effects of the boundaries are weak 
and the subchains in the blob follow excluded volume 
statistics. Namely, for each blob j = Cg\o, where gi is 
the number of monomers in a blob. The prefactor C is a 
constant, which depends on the temperature and the sol- 
vent. Thus, each blob contains gi = (%^)" monomers, 
and the number of blobs is rib = N/gi = N(^)i. Con- 
sidering the entropic elasticity of the polymer and the 
excluded volume repulsion between the blobs. The free 
energy of a linear chian confined in a nanochannel is given 
by Q 
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(3) 



where A and B are constants, the first term is the elas- 
tic free energy, and the second term represents excluded 
volume interactions. Minimizing the free energy with 
respect to R\\j leads to the extension j of polymer 
trapped in the nanochannel, we have 



Hi = (^) 1/3 D(N/ gi ) = {^ND^{Caf'\ 



(4) 



Correspondingly, the free energy at Rn ; is 



k B T 2 2 / 3 



A 1 / 3 B a ' 3 ND- 1 ' v (C*) 1 / v . (5) 



For a ring polymer confined in an infinitely long 
nanochannel of diameter D, there are always two folded 
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subchains inside the channel, which forms two string of 
blobs of size £b,r- In each blob, = Cg v r a with g r be- 
ing the number of monomers in a blob. The free energy 

due to excluded volume interactions is B^S- ^l 9 ^ . In 

H||, r 

contrast to the linear chain, the ring polymer can been 

considered as two subchains of length N/2, and thus the 

R 2 

elastic free energy is 2A ^ 2 ^'/2g ) • Therefore, the total 
free energy of a ring polymer confined in a nanochannel 
is 



ring 



knT 



= 2A 



m 



CL(N/2g r ) 



B^WltZ. (6) 

D 2 Rn r V ; 



The equilibrium size, -R||, r , at which T r i ng is minimized, 



= (fit WW 



Using £& >r = Cg^a and the Flory exponent v = 3/5, we 
further obtain 

R\Ur = (^) 1/3 ND^(Ca)^. (8) 

Interestingly, i£|L r is independent of the value o/£b. r . 

Obviously, Eqs. (j4]) and (jHJl derived from the Flory-likc 
theory are in good agreement with the Eq. @ from the 
scaling analysis. Moreover, Eqs. (JSj) and ((4]) also give the 
ratio 



as 



i?n 



Kl/3 



<4> 



0.630. 



(9) 



Compared with the value 0.718 for unconfincd chains, 
this ratio is smaller. 

However, the free energy of the chain at R^ r is a func- 
tion of r , which is 



k R T 



9_ = 3A y3 B 2/3 ND -4/3^/3-l/^ Ca y/u_ (1Q) 



It is difficult to determine £b. r . Taking into account ge- 
ometric average, 7r(£b ir /2) 2 = 0.57r(D/2) 2 , one obtains 

£,b,r = ^D. However, the geometric average here is 
questionable. Statistically, the blob should be spherical 
because the segments in the blob cannot feel the exis- 
tence of the wall. If we still use ^b.r = -^D, Eq. (fl"0)) is 
rewritten as 



k B T 2 2 / 3 - 1 /2^ 
As to the ratio the ratio 

rina 



A 1 l z B 2 l z ND~ 1 l v {Ca) 1 l v . (11) 
Eqs. (HTJ and © give 
= 1.782. (12) 
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FIG. 1: Longitudinal size Ru and the radius of gyration R g « 
along the channel axis of flexible self-avoiding linear chains 
as a function of aN(a/D)~~ 1 . Different polymer lengths N 
and channel diameters D are used. In the plots, we use v — 
0.588. The slopes, which stand for the prefactors of the scaling 
behaviors, are 1.367 and 0.399 for J?| and R g \\, respectively. 



For long chains partially confined in the nanochannel, 
it can escape from the channel by the entropic force. This 
pulling force can be estimated from the free energy and 
the longitudinal size, where prefactors of scaling behav- 
iors for the free energy and the longitudinal size are very 
important. The ratio of entropic forces for a ring polymer 
and a linear polymer is 



fr 



fi 



2 2/3+l/2„ = 2 ^/2 



i t tea r 



(13) 



where v = 3/5 is used. Even for one string of the long 
ring polymer, the pulling force acting on it is y/2 times 
as that for the long linear chains. This indicates that 
long ring polymer are easier to escape from the channel 
than long linear polymers due to smaller blob size in the 
channel for ring polymers. 

Fig. [1] shows the longitudinal size Ry and the radius 
of gyration R g \\j along the channel axis as a function 
of oN{ij / D)v~ l for flexible linear chains. Here, N is 
the chain length, a is the Kuhn length of the chain, and 
v = 0.588 is the Flory exponent in three dimensions 27 1. 
As expected, all the data for Ry and R g \\j collapse onto 
two different straight lines for different chain length TV 
and channel diameter D, respectively. This is in agree- 
ment with previous Monte Carlo simulation results [29(. 
The slopes, which stand for the prefactors of the scaling 
behaviors, are 1.367 ± 0.009 and 0.399 ± 0.003 for Ry 
and R g \\ t i, respectively. 

For flexible ring polymers, R\\ T an d R g \\,r as a function 
of <tN(<j/ D)"^ 1 for different chain length N and channel 
diameter D are plotted in Fig. O Interestingly, all the 
data for R\\, r and R g \\ t r also collapse onto two different 
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FIG. 2: Longitudinal size Ru and the radius of gyration R g ^ 
along the channel axis of flexible self-avoiding ring polymers 
as a function of aN(a/D)~ ! '~ 1 . Different polymer lengths TV 
and channel diameters D are used. In the plots, we use v = 
0.588. The slopes, which stand for the prefactors of the scaling 
behaviors, are 0.798 and 0.233 for 7?|| and R g \\, respectively. 



straight lines, respectively. The slopes are 0.798 ± 0.004 
and 0.233 ± 0.001 for R^ r and R g \\, r , respectively. 

These results demonstrate that a confined ring polymer 
behaves like a linear one, although they have different 
chain topology, in agreement with predictions in Eqs. ([2"]l. 
(HJ and ((8j) . Persson et al. [25| have measured the scaling 
of circular DNA extension with channel diameter. The 
scaling exponent they found is in agreement with our 
prediction in Eq. ((8|). 

The prefactor of the scaling behavior for the longi- 
tudinal size is a very important measure for the static 
and dynamics of the chain confined to the nanochan- 
nel. On the one hand, it is the extension factor of a 
chain in a nanochannel. Based on slopes of the curves 



[T] and [2j we have 



0.708 



0.584 and 



1.367 

2^32 ~ 0.584, which are quite close to the pre- 



in Figs. 

R oll,*- - 
Rglll 0.399 

diction in Eq. (jH]), reflecting that the ring polymer has 
less extension than the linear one with the same chain 
length due to the ring closure. At the same time, it 
demonstrates that the ring polymer with chain length 27V 
is 16.8% more extended than the linear chain of length 
N, and implies that for two overlapping polymers in a 
nanochannel each has a larger extension than that in the 
absence of the other. 
Recently, Levy et al. 



38| presented results concern- 



ing on partly folded DNA confined in straight nanochan- 
nels, where a dyed linear DNA was forced into the chan- 
nel in a folded configuration (where a portion near the 
middle of the molecule enters first) by applying an elec- 
tric field. They observed that the equilibrium extension 
per unit length for the folded portion was approximately 
30% larger than the equilibrium extension for unfolded 
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FIG. 3: Effect of the chain stiffness on 5"-^. Here, kb is the 

H n.i 

bending constant. The chain stiffness increases with kb- 



molecules. However, using confine spectroscopy, Persson 
et al. 



25| found that the extension at a given confine- 



ment for circular DNA is approximately 5% less than 
that for linear DNA for channel dimensions similar to 
those presented by Levy et al. [38[ . Obviously, there is a 
discrepancy between these two experiments. 

To solve this discrepancy, in the simulation we take 
into account the chain stiffness by by a harmonic angle 
potential Ub{&) = kb(9 — ir) 2 , where 9 is the bond angle 
between two consecutive bonds and kb is the bending con- 
stant. The chain stiffness increases with kb- Fig. [3] shows 



that 



first decreases from about 0.58 and then satu- 



rates to about 0.5 with increasing the chain stiffness. In 
above two experiments, different solvents are used which 
can lead to different stiffness of the DNA molecule. Thus, 
the discrepancy between these two experiments is from 
different chain stiffness. 

Nature not only imposes geometrical constraints on 
biopolymers by confinement through cell membrane, the 
cell nucleus or viral capsid, but also takes advantages 
of certain underlying chain topologies, such as the ring 
structure. In fact, E. coli has a rod-shaped geometry and 
its chromosome is not a linear polymer but a circular one. 
Based on Monte Carlo simulations, Jun and Mulder [l9[ 
have addressed a basic physical issue associated with bac- 
terial chromosome segregation in rod-shaped cell-like ge- 
ometry and found that two ring polymers segregate more 
readily than linear ones in confinement. According to our 
above theoretical analysis and simulation results, for ring 
polymers confined in a cylindrical nanochannel the blob 
size is smaller than that for linear polymers. This means 
that the pulling force for the ring polymer is larger than 
that for the linear one during the chromosome segrega- 
tion as shown in Eq. (|13|) , leading to a faster segregation 

m. 
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Conclusions. We investigate the chain conformation of 
ring polymers confined to a cylindrical nanochanncl using 
both theoretical analysis and three dimensional Langevin 
dynamics simulations. We predict that the longitudinal 
size of a ring polymer scales with the chain length and 
the diameter of the channel in the same manner as that 
for linear chains based on scaling analysis and Flory-type 
theory. Moreover, Flory-type theory also gives the ratio 
of the longitudinal sizes for a ring polymer and a linear 
chain with identical chain length. These theoretical pre- 
dictions are confirmed by numerical simulations. Finally, 
our simulation results show that this ratio first decreases 
and then saturates with increasing the chain stiffness, 
which has interpreted the discrepancy in experiments. 

Our results should enable a new understanding of the 
conformation statistics of confined ring biopolymers such 
as DNA. We believe that this work is important for 
understanding biological systems with more complexity, 
such as the behavior of DNA inside phages or the spatial 
organization of the bacterial nucleoid in E. coli. 
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